A flexible matrix Libor model with smiles 

Jose Da Fonseca* Alessandro Gnoatto^^ Martino Grasselli-'- 

March 1, 2013 



CN Abstract 

^^ We present a flexible approach for the valuation of interest rate derivatives based on 



^ 



Affine Processes. We extend the methodology proposed in Keller-Ressel et al. (2009 1 by 

C^ changing the choice of the state space. We provide semi-closed-form solutions for the 

^> pricing of caps and floors. We then show that it is possible to price swaptions in a multi- 

I factor setting with a good degree of analytical tractability. This is done via the Edgeworth 

CN expansion approach developed in Collin-Dufresne and Goldstein (2002|l. A numerical ex- 



,__i ercise illustrates the flexibihty of Wishart Libor model in describing the movements of the 

P^ implied volatility surface. 

d 

I 



> 

00 

^, 

rn 

o 



Keywords: Affine processes, Wishart process, Libor market model. Fast Fourier Transform, 
Caps, Floors, Swaptions. 



*Auckland University of Technology, Department of Finance, Private Bag 92006, 1 142 Auckland, New Zealand, 
email: jose.dafonseca@aut.ac.nz. 

^Universita degli Studi di Padova, Dipartimento di Matematica Pura ed Applicata, Via Trieste 63, Padova, Italy 
and Mathematisches Institut der LMU, Theresienstrasse 39, Munich Germany. E-mail: gnoatto@mathematik.uni- 
muenchen.de 

''Universita degli Studi di Padova, Dipartimento di Matematica Pura ed Applicata, Via Trieste 63, Padova, Italy. 
E-mail: grassell@math.unipd.it and Ecole Superieure dTngenieurs Leonard de Vinci, Departement Mathematiques 
et Ingenierie Financiere, 92916 Paris La Defense, France. 



JEL codes: G13,C51. 



1 Introduction 



In this paper we present a unified framework for the valuation of caps, floors and swaptions. 
These instruments are the most common derivative securities which are traded in a fixed income 



desk of a financial institution (see e.g. Brigo and Mercurio (2006 1). Practitioners usually price 



these products by relying on a Black-Scholes like formula, which was first presented in Black 



( 1976 1. The market convention of pricing caps and swaptions using the Black formula is based 
on an application of the Black and Scholes1 ( |1973 [) formula for stock options by assuming that 
the underlying interest rates are lognormally distributed. Remarkably, the use of this kind of 
formulae had no theoretical justification, since they involved a procedure in which the discount 
factor and the Libor rates were assumed to be independent so as to write the pricing formula as 
a product of a bond price and the expected payoff. The systematic use of this market practice 
ignited the interest of academics aiming at providing a coherent theoretical background. 



In a series of articles, [Miltersen et al.| ( |1997| l, [Brace et al] ( |1997| ), [Jamshidian] ( |1997| ) and 
Musiela and Rutkowski (1997J), provided these theoretical foundations, introducing the Libor 
and Swap Market Model. Following these papers a stream of contributions appeared, trying to 
extend the basic model to the case where the volatility of the underlying factor is stochastic. 



The most famous proposals on this side can be found e.g. in Andersen and Brotherton-Ratcliffe 
( [200T] ), [Wuand Zhang] ( [2006] ), [Joshi and Rebonato] ( p003| ), [Andersen and Andreasen] ( [2002] ), 



Piterbarg (2005a I, Piterbarg (2005b I. Other approaches explored different dynamics for the 
driving process with respect to the CEV or displaced-diffusion considered before for the Libor 



rate: for example Glasserman and Kou (20031, Eberlein and Ozkan (20051 introduced jump 
and more general Levy processes, allowing for discontinuous sample paths of the driving pro- 
cess. Another interesting approach is the one of IBrigo and Mercurio ( 2003| ) based on a mixture 
of lognormals. 



A typical problem in the previous approaches is that once the closed form solution for cap 
prices is found, to obtain an analogous result for swaptions it is customary to assume that the 
underlying (which is a coupon bond) behaves like a scalar process (typically again geomet- 
ric Brownian motion). This results in inconsistencies between the so-called Libor and Swap 
Market Models. Even more important, by assuming that the coupon bond is driven by a scalar 
process, we do not take into account the correlation effects among the different coupons, a key 
feature of a swaption which may be viewed also as a correlation product. This last remark is 



of paramount importance for practitioners (see e.g. the introduction of |Colhn-Dufresne and 



Goldstein (20021). 



In this paper we consider a new approach based on the stochastic discount factor methodol- 
ogy, where instead of modeling directly the Libor rate, one concentrates on quotients of traded 
assets (i.e. bonds). It has been first introduced by Constantinides (T992JI and then developed 



by Gourieroux and Sufana (201 1 1 in a spot interest rate framework and by Keller-Ressel et al. 



( 2009 1 in a Libor perspective. In this latter work they use affine processes on the state space 
M>Q as driving processes and provide a full characterization of the model, which allows them 
to provide closed form solutions for caps and swaptions up to Fourier integrals. This approach 
is very interesting and easily overcomes many difficulties which are to be faced in the compu- 
tation of Radon-Nikodym derivatives. 



We provide an extension of this approach, by considering affine processes on the state space 
Sj""*^, the set of positive definite symmetric matrices. This state space may seem awkward 
at first sight, but the processes belonging to this family admit a characterization in terms of 
ODE's which resembles the one found for standard affine models, an example being given by 
the famous puffie and Kan| ( [T9961 ) model. In fact, in |Cuchiero et al.| ( |2011| ) the authors extend 
to the set S^ (the set of positive semidefinite symmetric matrices) the classification of affine 

ffie et all t.003\ for the state space M^q x M""'^ introduced by 
( 1996 1. What is more, the state space S'^^ leads to stochastic factors which 



processes performed by 



Duffie and Kan 



are non trivially correlated. The most famous example of process defined in the set S'j'^ is the 
Wishart process, originally defined by Bru| ( |1991| ), introduced in finance by [Gourieroux and 



Sufana|p003| ) and then extensively applied in [Gourieroux and Sufanaj ( [20 10[ ), [Gourieroux and 
Sufana[ ( |20TT] ), [Da Fonseca et al.[ ( [2008] ), [Da Fonseca et al.[ ( [2007b| ), |DaFonseca et al] ( |20TT] ), 



Da Fonseca and Grasselli (201 1 1, among others. 



The interesting feature of our framework is the possibility to obtain semi-closed form solu- 
tions for the pricing of swaptions in a multifactor setting, which is a well known challenging 
problem. In fact the exercise probability involves a multi-dimensional inequality. There have 



been many approaches to simplify the problem: for example. Singleton and Umantsev (2002 1 
suggest an approximation of the exercise boundary with a linear function of the state variables. 



However, the most efficient approach seems to be the one of CoUin-Dufresne and Goldstein 



( 2002 1 which heavily uses the affine structure of the model and is based on the Edgeworth 



expansion for the characteristic function in terms of the cumulants. Since the cumulants decay 
very quickly the Edgeworth expansion for the exercise probability turns out to be very accurate 
and fast. 



The paper is organized as follows. In Section 2 we introduce our framework by recalling some 
useful definitions and results on affine processes. Section 3 investigates the case of the state 
space S^^ and presents the technical results. In Section 4 we focus on the pricing problem 
of the relevant derivatives. Caps and Floors are briefly treated since their pricing is now quite 
standard within the EFT methodology, while we devote more attention to the pricing of swap- 



tions by adopting the approach of CoUin-Dufresne and Goldstein (2002 1. Section 5 illustrates 
the flexibility of our framework through a numerical exercise. Section 6 concludes the paper, 
and we gather in the technical Appendices proofs and some remarks useful for implementation. 



2 Affine Processes on the set S^^ of strictly positive definite sym- 
metric matrices 

2.1 General results and notations 

To outline the setup we will consider affine processes taking values in the interior of the cone 
S'^. We will use the notations ipt{u) = 4'i't^u) and iptiu) = (t){t,u) so as to be consistent 



with Keller-Ressel et al. (2009 1. We will be employing a property of the functions defining the 
Laplace transform, that we report after the following 



Definition 1. ( Cuchiero et al. (2011 1, Definition 2.1} Let (il, J^, (J^t)t>o ' ^) ^^ a filtered prob- 
ability space with the filtration {J^t)t>o satisfying the usual assumptions. A Markov process 

S = (Si)^>Q with state space S^, transition probability pt{^o, A) = F{Ylt S A) for A G S'j", 
and transition semigroup (Pt)^>g acting on bounded functions f on Sj" is called affine process 

if: 

1. it is stochastically continuous, that is, linis^t ps{^o,-) = PtC^or) weakly on S^ Vt, x G 
Sj", and 



2. its Laplace transform has exponential-affine dependence on the initial state: 



Pte 



-Trf«S, 



OJ 



E 



-Tr[«Et 



-Fn 



e-T'-[<lp,(So,cie) = e- 



</.t{«)-Tr[Vt(«)So] 



Vt and So, u G S'^, for some function 



^>o 



x5j 



^>0 



and if) : 



^>o 



x5j 



(1) 



S' 



Having apphcations in mind, we will consider affine processes which are solvable in the sense 



of Grasselli and Tebaldi ( 2008p (who investigated affine processes on the more general sym- 
metric cone state space domain): this means that the state space that we will consider is the 
interior of S^, namely the cone of strictly positive definite symmetric matrices, denoted by 
5J^TM Solvability is important, in fact it ensures that the Riccati Ordinary Differential Equa- 
tion associated to the Laplace transform ([T]l through the usual Feynman-Kac argument has a 
regular globally integrable flow: this will be crucial to outline our methodology (see e.g. the 
proof of Theorem |4] in the sequel). 



The next property closes our survey on affine processes. It will be needed when we prove that 
the structure of the model is preserved under changes of measure. 



Lemma 2. nCuchiero et al. 



{2011 1 Lemma 3.2) Let S be an affine process on S't, then the 



functions (p and ip satisfy the following property: 



(f)t+siu) = (pt{u) + 4>s{i>t{u)), 

tpt+s{u) = 1ps{lpt{u)). 



'By analogy, the set of negative (resp. strictly negative) definite symmetric d x d matrices will be denoted by 
Sd (resp. Srf""). 



2.2 Examples 



The previous general framework may be quite abstract at a first sight, mostly because of the 
high technical level required to properly introduce the notion of admissibility and existence 
for affine processes, see Cuchiero et al. ( |2011 1. In this subsection we provide some examples 
which will illustrate some concrete applications. We start with the most important one, which 
will also constitute our main object of study in the numerical illustrations. 



2.2.1 The Wishart process 

We suppose that the process S is governed by the following (matrix) SDE: 



(2) 



which was first studied by |Bru ( 1991 1 and whose solution is known as Wishart process. We 
assume M, Q invertible and M negative definite so as to ensure stationarity of the process. 
Moreover we require ^Q^ = kQ^Q for a real parameter k > d + 1 to grant solvability (or 
equivalently to grant that DetCEt) > with probability 1). Under the solvability assumption 
Grasselli and Tebaldi ( 2008[ ) showed that the Riccati ODE corresponding to the characteristic 
function can be linearized and therefore admits a closed form solution. This is important in 
view of possible applications since in this case the functions (f) and tp in definition ([T]l are 
explicitly known: 

Proposition 3. Consider the process T, = (Si)Q^^^j, which solves the SDE ([2]). Then the 
conditional Laplace transform is given by: 



E 



-Tv\uT.i 



Ft 



-(/>^(u)-Tr[»/'r{M)St] 



where t := T — t. The functions (priu) cind V't(w) satisfy the following system ofODE's: 

dipr 



dT 
d(t)T 
dT 

which is solved by 



^r{u)M + M ' i)r{u) - 2^r{u)Q ' Q^r{u), V'oM = «, 



Tr 



kQ^ QiIjt-{u) , (/'ol'") = 



where 



and 



1pr{u) = {mpl2,T{u) + ^p22,T{u)) ^ (mpu^riu) + ^21,r(n)) 

M 2Q^Q \ 1 



1p21,Tiu) V'22,r(n) 



exp < r 



-M' 



4>t{u) 



-Tr 



log iuipl2,r{u) + ^22,t{u)) + M^ T 



(3) 

(4) 
(5) 

(6) 
(7) 
(8) 



Proof See GrasseUi and Tebaldi (2008 1. 



D 



The Wishart process constitutes the matrix analogue of the square root (Bessel) process. In fact 
we have that the matrix M can be thought of as a mean reversion parameter: this is evident 
from the Lyapunov equation defining the long-run matrix Soo, which is given by 



-kQ^Q = MSoo + Sooivr 



(9) 



The second way to appreciate the analogies w.r.t the square root process is to look at the dy- 
namics of the entries of the matrix process S. Concentrating on the main diagonal, in the 2x2 
case we have: 



dSn = {k {Ql^ + Ql^) + 2 (MnSn + M12S12)) dt 

+ 2aj^ {QudWl^ + Q2idWl^) + 2aj^ {QiidW2i + Q2idW^^) 
dT.22 = {1^ {QI2 + Q?2) + 2 (M21S12 + M22S22)) dt 



+ 2cjf {Qi2dWl^ + Q22dWP) + 2af {Qi2dW^^ + Q22dW^^) 



where we set 



o"ll cri2 
0"12 cr22 



S. 



(10) 

(11) 

(12) 



We refer to Da Fonseca et al. ( 2007a I for additional insights on the behavior of the Wishart 
process when aggregating its parameters. 

2.2.2 The pure jump OU process 

The procedure we adopt in this paper is general, meaning that we can consider different exam- 
ples of processes lying in the cone of positive definite matrices. In particular, we may consider 
the matrix subordinators proposed by [Bamdorff -Nielsen and Stelzer (2007 1, or jump-diffusions 
like in |Leippold and Trojani| ( |2010] ). In what follows we provide some examples with the cal- 
culations of the functions (j)r and ipr- 



Let us consider the SDE 



dSi = MSt + SiM ' + dLt 



(13) 



where M G GL{d) is assumed as usual to be negative definite so as to grant stationarity, and 
Lf is a pure jump process (compound Poisson Process) with constant intensity A and jump 
distribution v with support on S'^^ . The strong solution to this equation is given by: 

Et = e^*Soe*^^* + /* e^'^^'-^^dL.e^^^^*-''). (14) 

Jo 

We are interested in the computation of the Laplace transform of this family of processes: 



The functions and ip solve the following (matrix) ODE's: 



dr 
dr 



ls+\{o} ^ ' 



M'^t Mt 

e ue , 



The solution for the first ODE is given by: 

SO we can compute the Laplace transform by quadrature: 



-A/ \e 

/s+\{0} 



v{di). 



(16) 
(17) 

(18) 
(19) 



In the following we provide explicit computations by assuming some particular distribution 



u{-) for the jump size. The proofs of this formulae may be found in Gupta and Nagar (2000 1. 



For the sake of clarity, we specify that the Wishart distribution that we consider in the next 
sections are the classical distributions arising in the context of multivariate statistics. 



Wishart Distribution. Let J be the jump size. Consider the case J ~ Wisd {n, Q). Then 
we have 



0^(u) = -A / det (id + 2e*^^"ue*^"QJ ' ds + At. 
Non-Central Wishart Distribution. Let be J ~ Wisd (n, Q,^A), then we have 



exp < Tr 

+ Ar. 



^Q-^MM'^ + ^Q-^MM'^ Q-^'^e^^' "ue^^" + Q'^ 



Beta type I distribution. Let be J ~ /3^(a, b), then 

(priu) = -A / 1^1 (a; a + b; -e^^^ 'ue^^')ds + Ar. 
Jo 



(20) 



ds 



(21) 



(22) 



Beta type II distribution. Let be J ~ /3^^ (a, b), then 

Mu) = -A r ^^i^^{a■, -b+l{d+ 1); e*^"^ne*^^)d. + At, (23) 

Jo id(0j ^ 

where mFn, ^d{(^)y and ^(a; b; R) denote respectively the hypergeometric function of matrix 
argument, the multivariate Gamma function and the confluent hypergeometric function, see 



e.g. Gupta and Nagar (2000 1. 



3 A Libor model on 5^ 



To outline the general framework for Libor models, we start by considering a filtered measur- 
able space {Q,T,Tt) and a family of probability measures i^Tk)i<k<N- Under the measure 
Ftj^ we introduce a stochastic process S taking values on the cone state space S'j"''. At this 
stage the process may be a diffusion, a pure jump or a jump-diffusion process taking values on 
5j"+. Consider a discrete tenor structure = Tq < Ti < ... < T^ = T. We recall that the 
Libor rate is defined via quotients of bonds: 

1 fB{t,Tk-i) 



L{t,Tk):-- 



AT V B{t,Tk) 



1 



(24) 



where AT is assumed to be constant and AT = Tk — T^-i- The relation between the Libor 
rate and the forward price is given by: 



F{t,Tk-i,Tk) = l + ATL{t,Tk). 



(25) 



We proceed in full analogy with Keller-Ressel et al. ( 2009[ ) by extending their results to pro- 
cesses taking values on the cone of positive definite matrices. The intuition is simple: to build 
up a Libor model with positive rates, quotients of bonds should be strictly greater than one. On 



the other hand, a no-arbitrage argument (see e.g. Geman et al. ( 1995 1) implies that quotients 
of bonds must be martingales under the forward risk neutral measure indexed by the maturity 



of the denominator, so that the key ingredient in the approach of Keller-Ressel et al. (2009 1 
consists in the possibihty of constructing a family of martingales that stay greater than one up 
to a bounded time horizon. This will be possible thanks to the affine structure of the model, 
since in this framework bond prices are exponentially affine in the positive (definite) factors, as 
well as their quotients. 



3.1 Martingales strictly greater than one 

Let us first define the set 

Zt := {-u G 5d : E [e-Tii^s^]] < oo,VSo G S++] 
By the affine property of the process S we have 



E 



-Tr[«St 



-</.t(«)-Tr[t/-t(«)So] 



: [0, T] X Xr ^ K, 

V^:[0,T]xXT^5rf. 



(26) 



Within this setting we are able to construct martingales that stay greater than one up to a 
bounded time horizon T. 



Theorem 4, Let S be an affine process, and let u G It H S^ , then the process M" defined 
by 

M^ = exp {-(/>T-t(n) - Tr [V'T-t(«)Si]} (27) 

is a martingale and M" > 1 a.s. \/t G [0, T] . 

Proof. See Appendix. D 

Equipped with this positivity result, we can proceed by considering a tenor structure of non 
negative Libor rates L(0, T^) for k = {1, ..., N — 1}. Standard arbitrage arguments (see e.g. 



Geman et al. ( 1995 1) imply that discounted traded assets, in our case bonds, are martingales 



under the terminal martingale measure 



GA^(Pr^) VA;G{l,...,iV-l}, 



(28) 



B{*,Tn) 

where M^ (Ptjv ) denotes the set of martingales with respect to the forward risk neutral prob- 
ability Py^ . The idea in [Keller-Ressel et al.] (|2009[) is then to model quotients of bond prices 



using the martingales M'^ defined as follows: 

B{t,TN) 



M^' 



(29) 



B{t,TN^i) 



M] 



M]V-1 



(30) 



B{t,TN) 

Vt G [0, Ti] , ..., t G [0, Tn^i] respectively. As a consequence, the initial values of the martin 
gales Mq^ must satisfy the relation 

5(0, Tfc^ 



Mo"'= = exp {-<AtK) - Tr [V'tK)So]} - ^^^ ^^^ , 
for all k G {1, ..., N — 1}, so that it is possible to set uj^ = as we have M^^ = 1. 



(31) 



In the following proposition, we show that it is possible to fit (basically) any initial term struc- 
ture of bond rates. The state space we are considering offers a wide range of possibilities to 
perform this task. However, since we are interested in applications, we adopt the simplest 
choice directly coming from the scalar case and we focus on the particular (but realistic) case 
where all Libor rates are positive. 

Proposition 5. Let L(0, Ti), ..., L(0, T/v) be a tenor structure of positive initial Libor rates, 
and let S be an affine process on S^^. Define 

-Tr[MET] 



7s := sup E 



(32) 



lflT,> b(q't^ \ ^^^^ there exists a strictly increasing sequence of matrices (i.e. u^ -< lifc+i if 
and only ifuk — Uk+i G S^ ) ui <U2 < ■■■ < ujy-i -< in Xt Pi 5j~ and ujy = such that 



° ~ B{0,Tr,y 



VfeG {!,..., A^}. 



(33) 
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Conversely, let the bond prices be given by (|29p-(|J0|) and satisfy the initial condition {31 ). Then 
the Libor rates L(t,Tk) are positive a.s. Vt G [0, T^] and k G {1,...,A^ — 1}. 



Proof. See Appendix. 



D 



3.2 A fuUy-affine arbitrage-free model 

If we look at the definition of the Libor rate we realize that it is quite natural to require quotients 
of bonds to be driven by an exponentially affine function of the state: in fact, in this case also 
bond prices as well as forward prices will be affine functions. This is also in line with the 



previous approaches of Constantinides ( 1992 ) and Gourieroux and Sufana ( 2011 1 based on the 



stochastic discount factor. In other words, the approach of Keller-Ressel et al. (2009 1 is able to 
provide a fully affine structure^] To prove the affine structure or our model, we first show that 



under ([29|)-p0|), forward prices are of exponential-affine form under any forward measure. To 
do this, first we notice that in this framework quotients of bonds are exponentially affine in the 
state factors, so that also forward prices will be: for A; = 1, ..., A^ — 1 

Bit,Tk) _ B{t,Tk) B{t,TN) 






Bit, Tk+i) B{t, T/v) B{t, Tk+i] 

exp {Tr [(-V'TAr-t(nfc) + V'T^-t(ufc+i)) Si]} 

=: exjp{ATj^^tiuk,Uk+i) + Tr[BTj^^t{uk,Uk+i)^t\} ■ 



(34) 



With this result, we are able to show very easily that the model is arbitrage free, that is forward 



prices are martingales with respect to their corresponding forward measures (see Geman et al. 

Bi*,n 



( 1995 1) 



eM{¥Tj. 



(35) 



B{*,Tn) 

This comes from the fact that forward measures are related one another via the quotients of the 
martingales M": 



dF' 



Tk 



I .Ft 



F{t,Tk,Tk+i) B{0,Tk+i) 



M^" 



B{0,Tk) M^''+'' 



(36) 



\/k G {1, ..., N}. Then L{*, Tk) is a martingale under the forward measure Pt^+i since the 
successive densities from Pts,+i to ¥tj^ yield a telescoping product and a Py^^ martingale (see 



Keller-Ressel et al. (2009 1). More precisely: 

S(*,Tfe) 



l + ArL(*,Tfc) 



B{*,Tk+i) M 



— G M (Pt.^J 



smce 



Af-l 



M"* -r-r M"' 



l=k+l 



(37) 



(38) 



CoUin-Dufresne and Goldstein 



J2002| , who 



This is the reason why we will be able to apply the approach by ( 
originally started from an affine short rate in order to price swaptions: in fact, also in their framework bond prices 
are affine functions of the state variables. 



11 



Also, the density between the P^j. -forward measure and the terminal forward measure Ptjv is 
given by the martingale M"'= as indicated by ([29])-([30]): 



9Pt, 



\Tt 



BiO,TN) B{t,n) 



5(0, T,) * 



m; 



«fc 



Mo"'^' 



(39) 



9Pt^'-^' B{0,Tk) B{t,TN) 

In this arbitrage-free model with positive Libor rates, the affine structure is preserved: that is, it 
is possible to extend to the state space St^ the analogous result of 



Keller-Ressel et al. 



(20091. 



Proposition 6. Let the bond structure be defined through {29)-{30), where the process M"- is 
given by ^27\ . Then forward prices are exponentially ajfine in the state variable S under any 
forward measure. 



Proof. The result comes directly from formula (6.23) in Keller-Ressel et al. (20091 once the 
scalar product is replaced by the trace operator. D 



4 Pricing of Derivatives 

We now focus on the pricing problem for vanilla options like Caps, Floors and for exotic 
options like swaptions in the affine Libor model on S'^^ introduced in the previous section. We 
shall see that the pricing of Caps and Floors may be performed using standard Fourier pricing 



techniques as in Keller-Ressel et al. (2009), whereas, for the case of swaptions, we will resort 
to a quasi closed form solution. In fact, since the moments of the underlying affine process 
are known through its characteristic function, we can expand the exercise probability via an 
Edgeworth development, as shown in CoUin-Dufresne and Goldstein ( 2002[ ). This approach 
will lead to an efficient approximation that will avoid the numerical problems underlying the 



computation of the exercise probability in Keller-Ressel et al. (2009 1. 



4.1 Caps and Floors 

A Cap may be thought of as a portfolio of call options on the successive Libor rates, named 
Caplets, whereas Floors are portfolios of put options named floorlets. These options are usually 
settled in arrears, which means that the caplet with maturity T^ is settled at time T^+i. The 
tenor length AT is assumed to be constant. Since the two products are equivalent, we will 
focus on Caps. A Cap with unit notional has a payoff given by the following: 



AT {L{Tk,Tk) - KY 



k = l,...,iV-l 



(40) 



We rewrite the payoff of caplets as in [Keller-Ressel et aL] ( |2009[ ): 

AT {L{Tk, n) - K)+ = (1 + ATLin, n) - (1 + ATK))- 

with/C:=l + ArK. 



(41) 
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Thus we see that the caplet is equivalent to an option on the forward price. In order to avoid the 
computation of expectations involving a joint distribution, each single caplet is priced under 
the corresponding forward measure: 



:(Tfc,E:) = i?(o,r,,+i)E^^*+i 






/c 



icy 



with: 



Y := log 






ATN-n{uk,Uk+l) +Tr[BTr^_Tk{uk,Uk+l)^Tk] 



(42) 



(43) 



for ATj^_Tk{uk,Uk+i), BTf^^Tk{uk,Uk+i) defined as in (|34]|- The pricing problem can be 
solved via Fourier techniques through the |Carr and Madan (19991 methodology. Hence we 
have the following proposition, whose standard proof is omitted. 

Proposition 7. Let a > 0. The price of a caplet with strike K and maturity T^ is given by the 
formula: 

exp{— ac} 



CiTk,K) = B{0,Tk+i)- 



27r 



„—ivc 



=»(»'- (a+l)«)(^Tjv-Tj.(«fe,«fc+i)+Tr[BTjv-Tj.(«fc,Mfc+i)STj,]) 



(a + iv) {1 + a + iv) 



-dv, 
(44) 



where: 



c = log (1 + ATK) , 

ATN-Tk{Uk,Uk+l) = -(j)TN-Tk{Uk) +0Tjv-rfc(nfc+i), 
BTN-T^{uk,Uk+l) = -lJjTN-n(uk) + i^TN-n{uk+l)- 

In other words, pricing a Cap involves the computation of the moment generating function of 
e.g. the Wishart process, which can be efficiently performed through the linearization of the 
associated Riccati ODEs as explained in Proposition [3] The parameter a > represents the 
damping factor introduced by Carr and Madan ( |1999| ). We report in the Appendix B the explicit 
expression of the characteristic function involved in the pricing procedure. 



4.2 Swaptions 

The payoff of a receiver (resp. payer) swaption may be seen as a call (resp. put) on a coupon 
bond with strike price equal to one. We consider a receiver swaption that starts at Tj with 
maturity Tm, {i < m < N). The time-Tj value is given by: 



ET,iK,Ti,T„ 



Y^ ckB{Ti,n)-i 

\k=i+l / 



(45) 
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where 



Ck 



ATK ifi + l<k<m-l, 

1 + ATK if k = m. 



(46) 



Unfortunately, we face some difficulties if we try to adopt the Fourier technique that we em- 
ployed to price a caplet. To see this we look at the proof of Proposition 7.2. in [Keller-Ressel 
et al. (2009 1, which requires the computation of the Fourier transform of the payoff 



d(d+i) ^ 


V Cke'^^N-T^i^k,u^)+Tr[BT^.T,(uk,u,)T:T,] _ 


-1 


dvech{T^Ti)-, 




\k=i+l 


J 


(47) 



f{z) 



where for a symmetric matrix A, vech{A) stands for the vector in ]R^('^+i)/2 consisting in the 
columns of the upper-diagonal part of A including the diagonal. The problem is given by the 
presence of the positive part in the payoff function. To get rid of it, we should be able to find a 
value S such that 



k=i+l 



AT,^-TAuk,Ui)+Tr[BT,^-Tiiuk,Ui)f::\ 



1, 



(48) 



that is we should solve a single equation in d{d + l)/2 unknowns (the elements of S), which 
is highly non trivial when d > 1. Thus, pricing swaptions is challenging when we consider 



multiple factor affine models: this is a well known problem, see e.g. Jamshidian ( 1989 1 and 
Colhn-Dufresne and Goldstein|(|2002[l. [Keller-Ressel et al.] (|2009|l investigate the case d = 1, 



that is a Libor model driven by a (univariate) CIR process like in Jamshidian ( |1987| l. In that 



case, solving an equation similar to (48 1 is simple and the pricing of a swaption is only slightly 
more numerically complicated than the pricing of a Cap. As our purpose is to extend their 
methodology to a process with values in the set of strictly positive definite symmetric matrices 
we face a numerical difficulty related to the dimension of the state space. In order to solve 



this difficulty we follow CoUin-Dufresne and Goldstein (2002 1's methodology which strongly 
depends on the affine property of the process used to modelize the rates. As the processes we 
use have this affine property we can carry out the approximation for the swaption price pro- 
posed by these authors. Therefore, we can get around the dimensional difficulties posed by the 
process. 



We briefly recall the main results of CoUin-Dufresne and Goldstein ( 2002 1 to approximate the 
exercise probabilities for the swaption. We define the Tj-price of a coupon bond, fori<m< 
N, as follows: 

m 

CB{T,)= Y, CkB{Ti,Tk). (49) 

k=i+l 

Let us derive the general form of the pricing formula for a receiver swaption, for = To = t < 



u\±i, ±k} — B(Ti,Tf^) B(Ti,Ti) 



= exp {At^ -T,{Uk, Ui) +TT[BT^^Ti{uk,Ui)'ET,]} 
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8oiK,Ti,T„ 



E^ 



- I(f' rsds 



{CB{Ti] 



m 

J2 c,E^[e-/o'=^^^n(cB(T,)>i)' 



fc=i+l 

We switch to the forward measure as follows: 



So{K,Ti,TN)= Y. CkB{0,Tk)E^ 

k=i+l 



g- fo'^ rsds 

B{0,Tk) 



L(CB(T,)>1) 



B(0,T^)E^ 



5(0, T,) ' 



HCB(T0>1) 



fc=i+l 

- B{0,T,)E^^^ [l(^CBm)>i)] 

m 

= Y, ckB{o,n)¥Tj{CBm)>i)] 

k=i+l 

-B{0,T,)FtA{CB{T,)>1)]. 

The exercise probabilities PTfe [{CB{Ti) > 1)] andPr, [iCB{Ti) > 1)] do not admit in general 
a closed form expression but can be efficiently approximated thanks to an Edgeworth expan- 
sion procedure. Intuitively, the moments of the coupon bonds admit a simple closed-form 
expression in our affine framework, and these moments uniquely identify the cumulants of the 
distribution. One can expand the characteristic function in terms of the cumulants and compute 
the exercise probabilities by Fourier inversion. 



Using the notation of CoUin-Dufresne and Goldstein (2002 1 (their formula (5)) for the q — th 
power of a coupon bond we notice that, for i < ttt, < A'^: 

(CB(r,))» = (c,+iB(r,,r.+i) + ... + c„b(t.,t,„))' 



E 

jl,..;jg=l+l 



Cj, ■ ... ■ c,J X (B(r„r,J • ... • B{T„TjJ). (50) 



Now in our framework we have (see also formula (7.9) in Keller-Ressel et al. (2009 1) 



BiT^T, 






(51) 



for Z = 1, ...,q, meaning that we can rewrite the q — th power of the coupon-bond as follows: 



{CBm)r = Y 

iiv,i9=«+i 






m;- 






(52) 
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Recall, from pTll, that we have 



M'' = exp{-(/.Tjv-T,K) - Tr ['(/'t^_t,('UjJStJ} , 



for/ = 1, ...,q and 



m: 



T, 



exp{-(pTt,-TAui) - Tr [?/'T^r-Ti(^^)STJ} • 



In conclusion, the q — th moment under Ft^. has the following expression: 



^ (c,,.....c,JxE--. 

jl,...,jq=i+l 









E 



'-ji 



^J9 



) 



il,---Jq=«+l 



E'^^fe 



exp 



X]( -^Tn-T,' 



Uj 



Tr[V'T^_T,K)SrJ 



i=i 



+ ?('/'T^-T,(ni) + Tr [V'Tjv-T,(Ui)STj 



^ (cji • ... • cjj X exp <^ - ^ </'Tjv-r, (mj, ) + q^Tt,-% ( 



xE^^fe 



exp < Tr 



1=1 



(53) 



(54) 



(55) 



where the functions (p and tp are as usual the solutions of Riccati ODE's of the form (|4]), (|5). 
Once the first m moments under the corresponding forward measures are exactly determined, 
we can estimate the exercise probabilities Pt^ [{CB{Tq) > 1)] under each forward measure by 
relying on a cumulant expansion for Ft^ [CB(To)]. 

5 The Wishart Libor Model 



The aim of this section is to illustrate a specific choice for the driving process S. As in the 
general setup, we specify the process under the terminal probability measure Py^ . The example 



we choose is the Wishart process, which was already presented in section 2.2. 1 



dSi = {nn' + MEt + SiM ' )dt + VStdW/^Q + Q ' dW^ 



tTn I 



lT jTi/Tivl 



(56) 



Here W^ ^ denotes a matrix Brownian motion, i.e. a. d x d matrix of independent Brownian 
motions under the Pjv-forward probabihty measure. In the sequel we will write Wt for nota- 
tional simplicity. 
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In this section we show the impact of the relevant parameters on the impUed volatiUty surface 
generated by vanilla options for a Libor model driven by a Wishart process. With the aim 
to investigate some complex movements of the implied volatility surface, we first compute 
the covariation between the Libor rate and its volatility: this covariation is a crucial quantity 
allowing for the so called skew effect on the smile, in perfect analogy with the leverage effect 
for vanilla options in the equity market. 

5.1 The skew of vanilla options 

We want to compute the covariation between the Libor rate and its volatility, so we proceed 
to derive the dynamics of the Libor rate in the Wishart model. This may be done along the 
following steps: using the shorthand 

Bk ■■= BTf,^t{Uk, Uk+l) = -ll^Ti.-tiUk) + ^TN-t{Uk+l), (57) 

recall that we have: 



B[t,lk+i) 
In differential form, after dividing both sides by L{t, Tk, Tk+i) we have 
dL{t,Tk,Tk+i) 



(58) 



L{t,Tk,Tk+i) 

l + ATL{t,Tk,Tk+i) 



i[...]dt + Tr[BkdJ:t]). (59) 



L{t,Tk,Tk+i) 
To preserve analytical tractability, we freeze the coefficients and approximate as follows: 

I + ATL{t,Tk,Tk+i) _ l + ATL{0,Tk,Tk+i) ^. ^ 
L{t,Tk,Tk+i) L(0,Tfc,Tfc+i) 



Proposition 8. Under the assumption of frozen coefficients \60\ , the conditional infinitesimal 
correlation between the Libor rate and its volatility cannot be negative and is given by 



d {L{t, Tk, Tk+i),vol{L{t, Tk, Tk+i))) 

Tr [BkQ^QBkQ^QBkT;] dt 

Tr [QSfcS^TQT] ^Tr [^BkQ^QBkQ'^QBkQ^QBk] 



(61) 



Proof. See Appendix. D 

From the previous formula we realize that the matrix Q is responsible for the shape of the skew. 
We also have an indirect impact of the mean reversion speed matrix M coming from the term 
Bk which is the difference of two solutions of the Riccati ODE's Q and Q. The presence of 
S suggests that in the present framework the skew is stochastic. What is more, it can only have 
positive sign. 
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5.2 Numerical illustration with diagonal parameters 

The dynamics above show that the Wishart specification provides a very rich structure of the 
model. Since we want to get an understanding of the impact of different parameters we will 
look first at the case where all matrices are diagonal, which basically corresponds to a model 
driven by a two factor square root process (see e.g. Da Fonseca and Grasselli ( 201 1| )). 



We use the following set of parameters as a benchmark: 

_ / 3.75 \ / -0.3125 *1.0e~°°3 \ 

°~ y 3.45 y ' ~ y -0.5000 *1.0e-°°3 j ' 

^ / 0.034 \ 
\ 0.0420 J 

The impact of the Gindikin parameter k is quite easy to understand: the process acts by influ- 
encing the overall level of the surface. This is due to the fact that the higher k the lower the 
probability that the process S approaches 0. It is interesting to note that there is not only a level 
impact, but also a curvature effect, as we can see in Figure [T] 

[Insert Figure [T] here] 

Let us now look at the parameters along the diagonals of the matrices M and Q. The following 
claims may be easily checked by looking at the SDE's satisfied by the elements of S (see also 



Da Fonseca et al. (2007a I). Note that we assumed all eigenvalues of M to lie in the negative 



real line. 

— For I Mil I /^ (\) the surface is shifted downwards (upwards). 

— For IM22I /^ (\) the surface is shifted downwards (upwards). 

The impact is more evident for OTM caplets with short maturities. This is due to the fact that 
as the process decreases (in matrix sense) the probability that caplets with short maturities are 
exercised is lowered more than the analogous probability for longer term caplets. 

[Insert Figure |2]here] 

We then consider the impact of Qn, Q22- We have the following: 

— As Qii /^ (\) the surface is shifted upwards (downwards). In particular if we multiply 
Q22 by a constant c > 1, then the increment in the short term is higher for OTM than for 
ITM caplets. If c < 1 then the decrease is higher for short term OTM caplets, which is 
intuitive, given the discussion above. 

— The same impacts, with different magnitudes, is observed also for Q22- 

[Insert Figure |3]here] 



5.3 The term structure of ATM implied volatilities for caplets 

5.3.1 Diagonal parameters 

We proceed to consider the term structure of caplet implied volatilities. When the matrix Sq is 
diagonal, the impact of the elements of Q is the same: an increase in the absolute value of any 
element of Q will result in a steeper term structure of ATM caplet volatilities. 

[Insert Figure |4] here] 

Considering a model where Sq is a full matrix does not influence this result in a significant 
way. 

5.3.2 More complex adjustments: impact of off-diagonal elements 

To appreciate the flexibility of the Wishart framework, we focus now on the impact of the off- 
diagonal elements. We introduce off-diagonal elements in M and Q and look at the relative 
change in the short term smile (4 months) and the long term smile (32 months). We introduce a 
fully populated matrix Sq and look at the impact of M12 and M21 ■ Our experiments show that 
there is a symmetry between the sign of So,i2 (the initial value of S12) and M12, M21. More 
precisely, the implied volatility changes are as in Table [T] 



So, 12 > So,i2 < 



M12 > Increase Decrease 
M12 < Decrease Increase 

M21 > Increase Decrease 
M21 < Decrease Increase 

Table 1: Implied volatility changes: relation between So,i2 and M12, M21. 

The reason for this symmetry is to be looked for in the drift part of the dynamics of the single 
elements of the matrix process S. 

Next we look at the impact of Q12 and Q2i. To this end we model Q as a symmetric matrix 
and set Q21 = Q12 = P\fQuQr2 for a real parameter p. Also in this case we recognize two 
main shapes of the adjustment that we denote by S\^S2- 

We now proceed to perform other numerical tests which will show that our modeling frame- 
work has a certain degree of flexibility. For these tests we set: 

_ / -0.3125 *1.0e-003 q \ 

~ \^ -0.5000 *1.0e-°°3 J ' 

^ / 0.02 pVQUQ^A 

\ pvqi^q^ 0.0420 y 
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So, 12 > So,12 < 

p > Si 5*2 

p < 52 Si 

Table 2: Implied volatility changes: relation between 110,12 and p. 

so basically M is parametrized as before but Q is symmetric and equiped with a parameter p 
which summarizes the information on the off-diagonal elements. We require So = Sqo, where 
Soo is given by the solution of the Lyapunov equation Q. After that we perturbate Eq in order 
to include off-diagonal elements and set So,i2 = 5]o,2i = 2. We have a good degree of control 
on the term structure of ATM implied volatilities. In particular, we may have larger percentage 
shifts in the long-term w.r.t. the short-term ATM implied volatility, or, for p = —0.6 we may 
even reproduce a situation where the short term ATM implied volatility increases whereas the 
long-term ATM implied volatility decreases. 

[Insert Figure [5] here] 

If we adopt the same kind of parametrization for the matrix M by introducing a second param- 
eter p2, then we have further flexibility because we can impose many different combinations 
of p and p2. For example. Figure [6] shows that we are able to isolate an effect on the term 
structure of ATM implied volatility: in fact we have a moderate change for ITM caplets while 
OTM caplets are practically unchanged, but the shape of the term structure of ATM implied 
volatility is modified in a significant way. 

[Insert Figure [6]here] 

Finally, just for illustrative purposes we report a prototypical Caplet volatility surface generated 
by the model. 

[Insert Figure [7] here] 

As far as Swaptions are concerned an example of ATM implied volatility surface for different 
expiries and underlying swap lengths is given below. 

[Insert Figure [8]here] 

6 The Pure Jump Libor Model 

Finally, in this section, we would like to provide a second example for the driving process E, 
so as to let the reader appreciate the degree of generality of this framework. As in the general 
setup, we specify the process under the terminal probability measure Ptat- The example we 



choose is a matrix compound Poisson process, which was already presented in section 2.2.2 

dEt = MEi + EtM^ + JlJ"^^ . (62) 
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All assumptions presented in section 2.2.2 are in order. More specifically, we assume that L 



Pt^ 



is a compound Poisson process with constant intensity A and jump distribution taking values in 
5t^. As a specific example of jump distribution we choose the standard Wishart distribution. 



By recalling the results in section 2.2.2 we have that the solution for the function ipr{u) is 



1pr{u) 



e^'^^ue''^ 



whereas for (h-riu) we have 



(priu) 



-XI det{Id + 2e^^^'ue^^'Q] 'ds + Xr. 



(63) 



(64) 



In concrete pricing applications, the computation of the solution for (J)t{u) implies a numerical 
integration with respect to the time dimension. This numerical integration has an impact on 
the performance of the model which turns out to be slower than the Wishart Libor model. 
For illustrative purposes, we report an example for an implied volatility surface for caplets 
generated by the compound Poisson Libor model with central Wishart distributed jumps. The 
mean reversion matrix M and the jump intensity A are given by: 



M 



-0.0550 

-0.1760 



A = 0.1. 
As far as the jump size distribution is concerned, the parameters are the following: 

' 0.27 \ 



Q 



0.05 



n = 3.1, 



and the initial state of the process is 



_ / 1.875 0.6 
^~ \ 0.6 1.275 ^ 

[Insert Figure [9]here] 



7 Conclusions 



In this paper we presented an extension of the approach of Keller-Ressel et al. ( 2009 1 to the 



more general setting of affine processes on positive definite matrices. We showed that their 
methodology may be adapted to this state space in a straightforward way. What is more, it 
is possible to efficiently price European swaptions in this multi-factor setting by means of a 



cumulant expansion due to CoUin-Dufresne and Goldstein (2002). In doing so we are in front 
of a setting which is potentially able to capture correlation effects which can not be described 
by a single-factor framework. We provided numerical examples for the Wishart Libor model. 
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where the introduction of off-diagonal elements gives rise to new possibilities in the control of 
the shape of the implied volatility surface. 

Our contribution may be seen as a starting point for a description of market models in this state 
space, in consequence we believe that there are many possible directions for future research. An 
example is given by the problem of calibrating this family of models to real market data. As the 
structure of the products in the fixed-income market suggests, even in the plain vanilla case, we 
expect the objective function that should be minimized in the calibration procedure to be quite 



involved. Yet, some calibration results were obtained on equity derivatives in Da Fonseca and 



Grasselli (201 1 1 for Wishart based models so a calibration using interest rates derivatives might 
be feasible. Certainly, it will be a delicate issue and may constitute an interesting contribution 
by its own. Once the model is calibrated on vanillas, one could then further investigate the 
performance of the model on more exotic structures, like e.g. Bermudan swaptions and barrier 
options. Theses issues are left for future work. 
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Appendix A: proofs 

Proof of Theorem m 

For all li G Xt we have 

E [M^] = E 

and by the affine property we obtain 



-Tv[uT.t] 



< oo, 



E [M^\Tt] = E [exp {-(/>T-r(n) - Tr [iPt-t{u)^t]} |-^t] 
= E[exp{-Tr[nST]}|J't] 
= exp{-(j)T-t{u) - Tr [iljT-t{u)^t]} = M^, 

hence the process is a martingale. Now we show that M" > 1. Recall that by assumption 

It G Xt n S^^ and 

M," = E[exp{-Tr[nST]}|-Fi], 



so that if — Tr [uSy] > a.s. then we are done. We proceed as in Gourieroux and Sufana 



( |2003| l and apply the singular value decomposition to the negative definite matrix u, i.e. u may 
be written as: 

n 

T 



U 



1=1 



UjUi 



where Aj are the eigenvalues of u and Ui are the eigenvectors. By assumption S^ takes values 
in S'^^ , hence 



-Tr [uSt] = -Tr 



i=l 

i=l 
n 

y AjUj T^TUi > 



(65) 



i=l 



as we wanted. 



Proof of Proposition |5] 



We follow closely the proof in Keller-Ressel et al. ( 2009| ). By assumption, initial Libor rates 
are strictly positive, then 



i?(0,Ti) ^ B{^,T2) ^ ^ B{0,Tn) ^ ^ 



5(0, Tiv) B{0,Tn) 



B{0,Tn) 



Recall that we have 



E 



,-Tr[«iET] 



M^' = exp{-0T(ni) - Tr [^^1)^0]} 



B{0,Ti) 
BiO,TN)' 



(66) 



(67) 
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By the definition of 7x1 in (32 1, we have that if 72 = 00 then we are done, else we can claim 



that there exists an e > such that 7s — e > -dtWvH- Then we can find a matrix u s.t. 



E 



-TrfuSj 



BiO,TN)- 



> 7s - e > 



B{0,Ti 



B{0,Tn) 
In analogy with Keller-Ressel et al. ( |2009| ) we introduce the function 

/:[0,1] 



(68) 



^^E 



-Tr[gMEi 



(69) 



and we want to show that / is continuous. First, since S S Sj"*" and u G 5^~ we have that 
if M -< f then — Tr [uTjt] > — Tr [wSy], hence by monotone convergence we can conclude that 
/ is increasing. We now introduce an increasing sequence (an)neN /^ 1 and apply Fatou's 
lemma to obtain 



lim inf E 



-Tr[a,iMST] 



> E 



lim inf e"'^''!'^""^^! 



E 



-TrfuSr] 



implying that / is lower semi-continuous. Since / is also increasing we have that / is contin- 
uous. Now /(O) = 1 and /(I) > ^Ly^l , hence there exist some numbers = ^n < ^n-i < 
... < ^1 < 1 such that 

BiO,Tk) 



fiCk 



Mo«'=" 



VA;G {!,..., iV}. 



5(0, T^)' 

By setting u^ = ^fcU (for /c = 1, ..., A^ — 1) we obtain a sequence of matrices u^ -< ui^+i,Uk — 
Uk+i G 5^^ which allows us to fit the initial tenor structure of Libor rates as desired. Finally, 



we apply Proposition 1 and Lemma 3.2 (ii) in Cuchiero et al. (20 11 1 in order to obtain the last 
sentence of the Proposition [5j 

Proof of Proposition |8] 



In this section we proceed as in the proof of Proposition 4. 1 in Da Fonseca et al. ( 2008 1. Recall 



that Wt is a shorthand for W;^ . From (159|) it follows that 

dL{t,Tk,Tk+i, 
L{t,Tk,Tk+i) 



C\ i...)dt + 2JTr[QBkJ:BjQ 
:= C ( {...)dt + 2 JXr [QBk^BjQ-^jdWt] , (70) 




Tr 



QBkV^dWt 




where C was defined in ([60]) and the scalar noise driving the factor process may be derived as 
follows: 



dTr 



QBkT,tBkQ 



Tr 



QBk^Q ' QBkQ ' 



+ 2Tr 



QBkM^tBkQ' 



dt 



+ 2Tr 



QBk^/^tdWtQBkQ' 



{...)dt + 2 JTr [^tBkQ~^QBkQ^QBkQ^QBk 



:= {...)dt + 2 JTr [^tB^Q^ QB^Q^ QB^Q^ QBk]dZt. 



Tr [QBkQ^QBkV^dWt] 
VTr [^BkQ^QBkQ^QBkQ^QBk] 

(71) 
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The covariation between the noise of the Libor rate and its volatility is then given by 

Tr [QBkVWtdWt] Tr [QBkQ^QBkV^dWt] \ 

Tr [QBkJ^tBjQ-^] ' \/Tr l^tBkQ^QBkQ^QBkQ^QBk] / 

T,p,q,r,s QpqBgrV^rsdWsp) (Ea,6,c,d,e,/,g QabBbcQjdQdeBefVT.fgdWga 



(dWt,dZt 



/-^a,b,c,d,e,f,g,g,rBfeQedBcbQbaQaqBqrV'^rgV^gfdt 

Tr [QBk^BjQ-^] ^Tr [^BkQ^QBkQ^QBkQ^QBk] 

Tr [BkQ^QBkQ^QBk^t] dt 

Tr [QBk^tBjQ^] ^Tr [^tBkQ^QBkQ^QBkQ^QBk] 

Now we turn on the positivity of the skew. With the notation in the proof of Proposition [5j 
from ^k > ^k+i we have Uk ~< Uk+i and then Bk G S'^. In all terms in the numerator and 
the denominator we recognize congruent transformations of matrices in S'^ which leave the 
signs of the eigenvalues unchanged. The self-duality of 5j~ allows us to claim that all traces 
are positive, hence we are done. 



Appendix B: the characteristic function 

With the purpose of pricing caplets, we need to have a more explicit form for the characteristic 
function appearing in Proposition [7] Once we have this expression we can plug in the functions 
(pT {u) and t/jr (u) to obtain a closed form solution. The pricing problem will be then solved via 
FFT. Recall that we are considering the following expectation: 

gi(i,-(a+l)i)(Afc+Tr[BfcETj)" 



ip{v) = E" 



-k + l 



A-v-ia+mAk^vT^+i 



exp < Tr 



i{v- (a + l) i) Bk STfc 



(72) 



where 



Ak := -(l)TM-niuk) + (t>TM-n{uk+i); 

Bk ■■= -IpTN-Tkiuk) + IpTN-niuk+l). 



(73) 



As we computed the shape of the function (/)t{u) and tpriu) under the P^^^ -forward measure, 
we need to switch from the IPrfc+i to the P^^ -forward measure: 



^iiv-(a+l)i)Ak^FT^ 



Tk Ti[uT.Ti.] 



(74) 
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where the last equation follows from ([39]). Let us focus on the expectation which becomes: 



E^^iv 



exp I - (/>T^_r,K+i) - Tr [^T^_r,(nfc+i)STj 

+ (pTNiuk+i) + Tr [^Tjv(^fc+i)So] + Tr [uStJ | 

= exp I - (pTN-n (uk+i) + (t>TN (uk+i) + Tr [-0Tjv {uk+i)T,o] | 



pTr[(-VT^_Tjnfe+i)+«)STj 



exp 



I - (pTN-ni^k+i) + 0Tjv('"fc+i) +Tr[V'Tjv(^fc+i)So] 

<?^Tfe(-V'Tiv-TfeK+l)+uj -Tr VTfc(-V'TAr-Tfe(wA:+l)+njSo |. (75) 



Now, recalling the previous terms in front of the expectation in ( |74| ), we obtain the final expres- 
sion which is 



At: 



exp<i(z;- {a + l)i)\^-(l)Tj,^Tk{uk) + (pTN-niuk+i] 

- (t^TN-niuk+l) +(t)TN{Uk+l) 



Bh. 



-(pTk [ - 1pTM-niUk+l) + iiv - (a + l)i)( -tpTM-T^iuk) + VT^r-Tfc('"fc+l) 

+Tr[ipTMiuk+i)^o] 



Bu- 



Tr 



VTfc - ipTN-niuk+i) + i{v - (a + i)i) ( - VTjv-TfcK) + i'TN-niuk+i)) s, 
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Figure 1: Doubling k with respect to the basic case causes an upward shift of the surface. The 
plot represents the two smiles (4 months and 32 months) for the basic {k = 3) and the modified 
case (k = 6). 

Comparison w.r.t. the basic case M, =-0.0005625 
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Figure 2: Impact of Mu. Mn is negative and the present image shows the effects on the two 
smiles (4 months and 32 months) we obtain when we multiply it by a constant c = 1.8 
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Figure 3: Impact of Qu. Qu is positive and the present image shows the effects on the two 
smiles (4 months and 32 months) we obtain when we multiply it by a constant c = 2 
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Figure 4: Impact of Q on the term structures of ATM implied volatilities. Here we consider 
Qii and multiply its value by a constant c = 1, 1.5, 2 so as to get the values in the legend. 
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Figure 5: The images above highlight the flexibility of the Wishart Libor model. We are able 
to impose different patterns to the term structure of ATM implied volatility. On the top we 
have the smiles and on the bottom we observe the relative changes of the smiles, i.e. for every 
point of the smiles we calculate the quantity (cj}7„^^; - aZfuai) /'^Zuiai- Notice in particular 
the situation on the left side, where we observe around 5% (ATM) an increase of the short term 
smile and a decrease on the long term. 
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Figure 6: Impact on the implied volatility surface when both M and Q are parametrized as 
symmetric matrices. Notice the level around 5%, corresponding to ATM. This shows that if 
we parametrize both M and Q via p, p2 we have a flexible setting which is controlled just by 
two parameters that allow us to perform different combinations. In particular p and p2 have 
opposite impacts in the present example (/) > whereas p2 < 0), meaning that we have a good 
degree of control. 
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Figure 7: Caplet Implied Volatility Surface generated by the Wishart Libor model 
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Figure 8: ATM Swaption Implied Volatility Surface generated by the Wishart Libor model 
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Figure 9: Caplet Implied Volatility Surface generated by the compound Poisson Libor model 
with Wishart distributed jumps 
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